A nearly closed ballistic billiard with random boundary transmission 
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A variety of mesoscopic systems can be represented as billiard with a random coupling to the 
exterior at the boundary. Examples include quantum dots with multiple leads, quantum corrals with 
different kinds of atoms forming the boundary, and optical cavities with random surface refractive 
index. We study an electronic billiard with no internal impurities weakly coupled to the exterior by 
a large number of leads. We construct a supersymmetric nonlinear <r-model by averaging over the 
random coupling strengths between bound states and channels. The resulting theory can be used 
to evaluate the statistical properties of any physically measurable quantity. As an illustration, we 
present results for the local density of states. 



When the size L of a two-dimensional (2D) mesoscopic 
structure is reduced beyond the elastic mean free path 
I of the electrons, all the scattering takes place at the 
boundary. An electron may escape through the bound- 
ary, get reflected back specularly, or participate in both 
processes. Besides quantum dots hi, which can exhibit 
both diffusive and ballistic behavior|2>l3> there are meso- 
scopic billiards which are ballistic by construction. Ex- 
amples include quantum corrals (QC) 0], optical corrals 
0, optical resonant cavities 0, and the artificial atoms 
proposed in Ref. 0- 

A quantum dot has boundary losses, unless the confin- 
ing potential is chosen to be infinitely high. Such losses 
are usually modeled by coupling to a number (possibly 
infinite) of open channels jg, although the precise de- 
tails of coupling are generally unknown. In this paper 
we focus on internal ballistic dynamics of a clean circular 
dot which is nearly closed, i.e., it is weakly coupled to 
a large number of leads. The crucial ingredient in our 
model is the randomness of the transmission coefficients. 
It enables us to carry out an ensemble average using the 
supersymmetry method @], and, as a side-benefit, acts as 
a natural regularizer, helping us avoid the technical diffi- 
culties of previous supersymmetric approaches to closed 
ballistic systems IsJLLlJ]- The resulting theory is a "sur- 
face" cr-model, which resembles the conventional diffusive 
tr-model 0, but has the "diffusion" modes confined at the 
boundary of the dot. These can be associated with clas- 
sical whispering gallery trajectories, which run along the 
walls of the system and are known to str ong ly influence 
transport through mesoscopic structures [jllj. 

Furthermore, as pointed out in Ref. ^3], the ballistic 
analogs of density relaxation modes in diffusive samples 
originate from trajectories which remain close to each 
other in configuration space. Given the almost closed 
nature of our dot, the whispering gallery modes (WGM) 
fill are expected to impact the long-time characteristics 
of the internal dynamics and to play a major role in quan- 
tum interference effects. By contrast, the "star-like" and 
"asterisk" trajectories, which approach the lead mouth at 
small angles to the boundary normal and are more likely 
to exit the billiard, cannot contribute 0] to the response 



functions. In other examples 0, Q] the electrons do not 
decay into well-defined leads, but because of the signif- 
icant similarities between the three nanostructures our 
approach will apply, mutatis mutandis, to them as well. 
To be specific, the long-lived modes which play a chief 
role in their internal dynamics have the same classical 
origin. 

A quantum dot in which electrons are confined by a 
hard wall potential but can escape into leads can be 
described by the non-Hermitian effective Hamiltonian 
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Here, Hq is the Hamiltonian of the closed quantum dot, 
and 5c n is a surface <5-function with uniform support on 
the crossection C n which separates the nth lead from the 
billiard. It is important to note that the leads have been 
"integrated out". In other words, one replaces an open 
dot with plane wave boundary conditions in the asymp- 
totic region and Dirichlet conditions along the rest of the 
boundary, with a closed dot and simplified boundary con- 
ditions The operators B n are defined through their 
action on an arbitrary function 
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where 7„ < 1 is a coupling coefficient, Vn' is a trans- 

(i) 

verse velocity, and ip,i is a normalized eigenfunction of 
transverse motion for zth channel in nth lead. 

Thus, the solution of the original open dot problem 
reduces to the solution of 

(E-H eff )G R < A (ry)=d(r-r>), (2) 

for the retarded (advanced) Green's function G R ' A , with 
Neumann boundary conditions (x„v r ) G R A (r,r')| c = 
0, where x„ is a unit vector parallel to the waveguide 
walls, and v r is the velocity operator. It acts from the 
side of the lead. 

Just as in studies of chaotic scattering 0], the key tool 
in calculations of response functions of ballistic (or disor- 
dered) systems with leads is the operator (E — H e /f)~ . 
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In the framework of the supersymmetry method it can be 
treated in the same way as (E — Ho) 1 , as was done in 
Ref. 0| > where nonperturbative calculations were car- 
ried out for disordered dots. Note, that we pursue the 
non-universal, i.e. non-"random matrix theory" regime. 

For a ballistic dot the dynamics is governed by the 
operator: Ho = (p — (e/c)A) 2 , where we assumed no 
potential and a constant magnetic field. The latter is 
introduced to break time reversal symmetry, which sim- 
plifies the application of the supersymmetry method, and 
makes the illustration of our approach more transparent. 

We intend to calculate averages over the Gaussian 

(i) 

distribution of dimensionless coupling coefficients 7„ . 
The relation of these coefficients to sticking probabili- 
ties, transmission coefficients and other commonly used 
parameters can be found in Ref. in the context of 
the Hamiltonian approach to chaotic scattering. More 
insight on the physical meaning of these coefficients can 
be gained from Refs. [f^ |. 

Before commencing our supersymmetric derivation, we 
make a few simplifying assumptions. These assumptions 
can all be relaxed without affecting the physics, but are 
needed to simplify the technical details. We allow only 
one open channel in each lead and express the coupling 
coefficients as a sum of constant and stochastic parts: 
7n = 7 + In- For the statistics of 7„ we assume that 
(7n) = 0, (7n7m) = x 2 S nm , and that all higher moments 
factorize into second moments. Next, we "eliminate the 
leads" 0] , passing to the Hamiltonian given by Eq. Q . 

Upon introduction of four-component supervectors 
* (r) T = {Si (r) , xi (r) , S 2 (r) , Xi ( r )}, and supermatri- 
ces L = diag {f , f — 1,1}, A = diag {1, 1 — 1,-1} (see, 
for example, Ref . 1 171 1. the supersymmetric generating 
functional reads |l£| 



(Z [J])~ n = J dV*dH>e- c W <V C * [ * ] )_ , (3) 

where (...)= indicates averaging over random couplings 
to the leads 

£[#]=« J ¥(r) HL^{v)dr 
~ N f 
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with H = H0I4 + ieA, (e is infinitesimally small) Ho — 
— V 2 /2m — E, cp n (y) = y/2/d n sin (ny /d n ) (for hard-wall 
lead of width d n ). Here J c stands for a double integra- 
tion over y n and y' n , the transverse coordinates along the 
crossection C„ (perpendicular to the walls of a waveg- 
uide); the product dy n dy' n will be omitted in what fol- 
lows. 



Averaging over 7„ produces: 



Then, using the Hubbard-Stratonovich transformation, 
we decouple the "interaction terms", introducing super- 
symmetric fields Q n (yn,y' n ) defined at each crossection 
C n . To simplify the form of the action we further assume 
that the leads have identical width d and are attached ev- 
erywhere along the perimeter of the dot. We choose the 
total number N = 2irR/d (R is the radius of the cir- 
cle enclosing the dot) of leads to be large, and therefore 
d <C R. We also assume identical transverse velocities: 
v n = v. In the limit N — + 00, the density field Q becomes 
continuous. It turns into a function of a single variable - 
the polar angle 9. 

After the Gaussian integration over \& (r)-variables, an 
important supermatrix which needs to be determined is 
the effective Green's function Q (r, r'). It satisfies 
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Q(9))g(r,r')=iS(r-r'). 
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where Q (9) = Q (9) - (j/2xmd) A and a 
generating function takes the form 

W]> 7 = / 'DQe F M +F 'W, 
with the free energy 



F [Q] = Str J drdr' 
5(r-R) 
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*(r-r / ) + ln-»a- 1 (r,r')L (4) 



and the symmetry breaking term F e [Q] = 
-Str/ln(I 4 + eA<?(r, r'))drdr'. The functional in- 
tegration over Q (9) is performed in saddle point 
approximation, which requires the solution of 

TA1 

Q sp (9) = —g(R,R,9,9,Q sp (9)). (5) 

Assuming the solution to be diagonal and coordinate in- 
dependent, we arrive at usual structure of the saddle 
point: Q sp (9) = QoA. In order to analyze the fluctua- 
tions, the constant Qo and the diagonal Green's function 
supermatrix Q sp (r, r') are necessary. Thus, we mapped 
the original problem with random boundary condition 
onto the effective problem specified by the differential 
equation: 



(V 2 - S 2 )^ p (r,r',9,9') 



2miS (r - r') 8 (0 - $') 



(6) 



and the uniform boundary condition 



^G sp ( r ,r\9,9')\ s ^ i b -^g sp (r,r',0,0')l S - , (7) 
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where s 2 — —2mE, b = m 2 avx, and S~ is the inner 
surface of the dot. For both inner (r < R) and outer 
(r > R) domains the solutions are readily obtained and 
matched together. Below we will only need the solution 
for the inner domain, which can be written in terms of 
the modified Bessel functions /„ and K n of nth order 

G sp ( r ,rV ) = ^£ Msr<) 
l 

x{ a; / i (5r>) + ^( S r > )}e li ( 9 - 9 '), ( 8 ) 

and coefficients a/ are given via 

-ibQ' Q K t (sR) + sRK[ [sR) 



at 



ibQ' Ii (sR) - sRI{ {sR) ' 



Then, Qq is obtained from the stationary point condition 
(Eq. (0). Dropping the imaginary part of Qq, since it 
can be absorbed into E, we get 



E 



2tt (l 
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(9) 



Next, we set g = sR, f = bQo, and evaluate the sum over 
n in Eq. © asymptotically in the following limit: j> 1, 
f/g ~ 1. We replace the sum with the integral, switch to 
the new variable fi = v/ 'g, u se uniform expansion for the 
Bessel functio n !„ [20] and expand the integrand in 

l/g; see Ref. ^| f° r the details. After the substitution 
g — > — ig (s —* —ik) the left-hand side of Eq. © becomes 
f 2 — g 2 to the leading order in l/g. Therefore, we 
have an algebraic equation for / (or Qo), which can be 
solved with the help of the condition k 2 — (ir/d) = m 2 v 2 

At this point we introduce several relevant energy 
scales in "natural" units h = c = 1: Thouless energy 
Et — kj (mR), mean level spacing of a closed dot A = 
1/ (mR 2 ), and a total resonance width T — 7/ (mdR) 
related to the modal decay rate of a dot having identi- 
cal coupling coefficients 7 to all N channels. In our al- 
most closed dot we have the following hierarchy of scales: 
Et 3> A 3> T. Therefore, the set of dimensionless 
parameters specifying the problem completely is given 
by g = Et/A = kR 3> 1 (dimensionless conductance), 
M. = r/A <C 1 (modal overlap) together with 7, x and 
v. Any other constants which enters Eqs. QEJ, can be 
expressed in terms of these five. Hereafter we continue 
to use the old set of parameters, including / and Qq, to 
keep the notation compact. The energy scales we just 
specified are to be used in comparison of our predictions 
to experimental and numerical results. 

Now, we turn to the fluctuations around the saddle 
point, which can be decomposed into a transverse piece 
SQ^ (along the saddle point manifold and a longi- 
tudinal piece 8Q( l > (orthogonal to the saddle point man- 
ifold). The part of the action corresponding to the SQ^' 



(anticommuting with Q) is given by [2l| 
F t [SQ] = -m 2 a j (0)) ' d0 + ( 
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Q (Q) Q (-Q) 5Q (t) (9) &Q (t) (9') d9d9'. 



Expanding in angular harmonics SQ^ (9) = 
ex P {^9} /2-7T and using the Ward iden- 
tity (relation between product and difference of Q (Q) 
and Q (— Q)) it is possible to show that the massive term 
in Ft[SQ] is proportional to mRj/TrxQo <g; 1. Then, 
applying the same technique, as in the solution of Eq. 
(EJ, we find that, to the leading order in l/g, the free 
energy is quadratic in SQ^ (9) for vanishingly small e, 
and arrive at 



-Do 



mxvan 



g(2g 2 + .f 2 ) 
4Qo J P^l 2 ^' 



To finish the construction of the nonlinear cr-model, we 
integrate out the longitudinal modes, which decouple 
at this order in l/g from the transverse ones, and set 
«JQ« = Q. Next, we expand the symmetry breaking 
terms F e [Q (Q)] to the lowest order in e. The result is 
given by (Z [J]} = j DQe~ F ^\ with the free energy 



F\ 



Str / drdr' i D 



dQ 
09 



8{r-R) 



5(r-r') 



d9"eAQ (9") a (r, R, 9";r', R, 9") 



(10) 



a (r, R, 9"; r', R, 9") = i^Q sp (R, 9", r') G sp (r, R, 9") 

where r> (r<) is a maximum (minimum) of |r| and 
|r'|. The supermatrix Q satisfies a nonlinear constraint 
S 2 = Qoli an d can be parametrized as suggested in Refs. 
0, E| for the diffusive case. The n-point correlations 
can be generated from the functional given by Eq. (11(11 . 
which is the main result of this paper. Just as in case 
of the supersymmetric nonlinear cr-model of Ref. , the 
diffusion modes clearly play an important role in the su- 
perintegrals representing correlators. 

A physically measurable quantity which does not 
depend on fluctuations around the saddle point is 
average local density of states (LDOS) (p(r)) 7 = 
— (1/n) 3 (G (r, r', E)) . Our Indeed, this one-point 
function neither requires the knowledge of the A- 
like structure of the saddle point manifold, nor its 
explicit parametrization. It can be shown that, 
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3(G(r,r',£)). 



ic 



)>c 



where (• ■ ■ }q stands for integration with weight 
exp{— -F [Q]} (Eq. l(Toll ). This integration reduces 
to evaluating the integrand at the saddle point. Most 
conveniently, the average LDOS can be calculated 



via regularized resolvent K (r, s 
g (r,r,e,9,s 2 ) 0: 7r<p(r)) 7 = 
The result reads 



s— — ik 



-J 



E 



In («#) 



f 2 i?Ag) + g 2 (iMY 



(ii) 



where K = r/R. This expression for *zsK (s 2 ) can be 
further simplified for the limiting cases: (i) k -C 1/g, 
when the main contribution to sum in Eq. 11 111 comes 
from the terms with small n, and (ii) k > 1/g, when we 
can employ the uniform expansions of /„. The results 
are summarized in Fig. 
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Figure 1: LDOS plotted as a function of k in units of m/ix 
for two sets of parameters: g = 50, M — 1/12 (dashed line); 
and g = 200, M = 1/3 (solid line). In both cases the rest of 
the parameters were fixed x = 0.1; 7 = 0.2; v = 100; 

Thus, we have constructed a nonperturbative theoret- 
ical framework to analyze one particular realization of a 
whole class of nanostructures: a nearly closed system 
with ballistic internal dynamics interacting randomly 
with the outside world through the boundary. Our ap- 
proach introduces a natural regularizer, which enables us 
to circumvent the conceptual difficulties of previous ap- 
proaches to closed ballistic systems 0] . We find that the 
resulting theory, encapsulated by Eq. ltlfl|l . can be char- 
acterized by diffusive modes confined to the boundary 
and interacting nonlocally with the interior (see the last 
term of Eq. IjlOJl V The supersymmetric functional was 
constructed with the help of large angular momentum 
modes identified as the WGM. These modes are expo- 
nentially less likely to escape 0] compared to the modes 
with incidence directions close to the lead normals, and 
consequently, they dominate response functions at large 



times. Our framework should allow us to compute the 
statistical properties of any physically measurable quan- 
tity, though technical difficulties may impose strong lim- 
its. It should be clear that our approach is also applica- 
ble, with minor modifications, to other examples belong- 
ing to this class of systems. 

Finally, the extension of this approach to generic bil- 
liards with smooth walls (to be published elsewhere) is 
also possible, although it is more technically involved 
[jfjj . At first sight the non-linear supersymmetric a- 
model (NLSerM) for the rough billiards proposed in Ref. 
[2^ | looks very similar to the NLScMs we derived here for 
open circular billiard (and for open rough billiard in Ref. 
|2fij). However, there are several differences between two 
models, which can be summarized as follows. The dif- 
fusion and (one-dimensional) localization in angular mo- 
mentum I space described in Ref. |22j is guaranteed by 
small changes in I as the particle bounces off the walls. In 
our case, because of the sharp edges of the region which 
connects the leads to the dot the WGM trajectories may 
have much larger I increments along the way. As a result 
our model describes diffusion and localization in position 
(angle 9) space rather than angular momentum space. 

Another issue is the role of electronic interactions. One 
of the possible ways to take them into account in diffu- 
sive and ballistic systems with large dimensionless con- 
ductance, is to use a "Universal Hamiltonian" [3], which 
was shown to be the renormalization group fixed point for 
weak interactions j3 EB| • We hope to extend our anal- 
ysis to the interacting ballistic case by using the large- TV 
approach of Ref. j2^| . We leave these questions for future 
work. 
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